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FOREWORD 


This  report  was  prepared  by  the  M.I.T.  Gas  Turbine  Laboratory 
under  Contract  No.  N000LA-76-C-0253  and  covers  the  third  year  of 
activities  under  this  contract. 

The  principal  investigator  is  Professor  Jean  F.  Louis.  The 
contract  monitor  is  Mr.  James  R.  Patton,  Jr.,  Power  Program,  Office 
of  Naval  Research,  Department  of  the  Navy. 


ABSTRACT 


In  the  third  year  of  the  contract,  further  advances  were  made  towards 
the  goal  of  gathering  the  heat  transfer  and  aerodynamics  flow  data  necessary 
for  a good  understanding  of  the  performance  of  film-cooled,  highly-loaded, 
transonic  turbine  blading. 

The  MIT  cascade  blowdown  facility  now  fully  operational  was  used  in 
evaluating  the  heat  transfer  performance  of  the  four  blade  profiles 
designed  in  the  first  year  of  the  program.  The  results  show  that  the  level 
of  turbulence  is  an  important  parameter  in  determining  heat  transfer  in 
transonic  cascades.  It  also  shows  that  the  heat  transfer  to  the  trailing 
edge  of  the  blades  is  very  high  being  about  75%  of  the  heat  transfer  to 
the  leading  edge.  A comparison  of  the  Nusselt  number  calculated  from 
heat  transfer  measurements  with  the  Nusselt  number  obtained  by  a prediction 
method  using  the  pressure  distribution  shows  good  correspondence.  The 
variation  of  average  Stanton  number  over  a range  of  Mach  numbers  shows 
that  the  reference  blade  has  the  most  superior  heat  transfer  performance. 
Preliminary  data  has  been  obtained  on  the  off-design  performance  of  the 
blades  and  full  scale  tests  are  underway.  Comparative  studies  show  that 
about  21%  less  heat  needs  to  be  taken  out  by  internal  cooling  if  one  stage 
of  a transonic  turbine  is  used  to  replace  two  moderately  loaded  subsonic 
stages  which  produce  the  same  output,  have  the  same  inlet  stagnation  con- 
ditions, have  the  same  mass  flow  and  the  same  tip  speed. 

This  demonstrates  one  of  the  potential  advantages  of  transonic 
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STUDIES  ON  TRANSONIC  TURBINES  WITH 


FILM-COOLED  BLADES 
A THIRD  ANNUAL  REPORT 

H.  0.  Demuren,  N.  Adams,  F.  Hajjar,  0.  M.  Amana,  J.  F.  Louis 

Introduction 


The  increasing  interest  in  Transonic  turbines  shown  over  the  last 
few  years  is  a reflection  of  the  need  to  increase  the  work  output  obtained 
per  unit  of  cooled-blade  surface.  Much  of  the  increase  in  turbine  efficiency 
obtained  over  the  last  two  decades  has  been  due  to  our  ability  to  increase 
the  turbine  inlet  temperature  through  the  use  of  more  efficient  cooling 
techniques.  As  more  and  more  coolant  mass  flow  is  obtained  from  the  com- 
pressor, however  , a point  is  reached  when  the  losses  introduced  by  the  use 
of  coolant  mass  flow  more  than  offset  the  increase  in  efficiency  attainable 
by  an  increase  in  the  turbine  inlet  temperature.  At  this  point,  the  need 
arises  to  optimize  the  work  output  per  unit  of  cooled-blade  surface.  This 
requires  the  use  of  transonic  blading. 

The  goal  of  these  studies  has  been  to  gather  the  necessary  heat 
transfer  and  aerodynamic  data  which  would  be  of  use  to  the  designer  of 
high  performance,  film-cooled,  highly-loaded  transonic  turbine  blading. 

In  the  first  and  second  years  of  the  contract,  the  following  key  tasks 
were  performed: 

(a)  The  gathering  of  the  film  cooling  effectiveness  data  in  the 
transonic  range  1.5<M<1.4  using  the  shock  tunnel  for  both 
slot  and  holes  injection. 

(b)  The  gathering  of  heat  transfer  and  aerodynamics  data  at  the 
trailing  edge  of  a transonic  blade, and  the  use  of  this  data  to 
formulate  a model  for  the  heat  transfer  and  aerodynamic  flow 
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at  the  trailing  edge  of  transonic  blades. 

(c)  The  design,  manufacture  and  commissioning  of  the  MIT  cascade 
blowdown  facility. 

(d)  The  design  and  aerodynamic  testing  of  four  transonic  blade 
profiles  in  a conventional  wind  tunnel  at  VKl. 

(e)  An  estimation  of  the  coolant  flow  requirements  of  a transonic 
turbine  and  a comparison  between  the  coolant  flow  requirements 
of  the  transonic  turbine  with  that  of  a subsonic  turbine  of 
equal  work  output. 

In  the  third  year,  the  following  tasks  were  accomplished: 

(1)  Aerodynamic  and  heat  transfer  tests  were  performed  for  the  four 
transonic  blade  profiles  designed  in  the  first  contract  year, 
at  zero  angle  of  incidence  in  the  blowdown  facility. 

(2)  The  instrumentation  in  the  blowdown  facility  was  updated. 

(3)  A preliminary  evaluation  of  the  off-design  performance  of 
the  reference  blade  profile  has  been  carried  out  and  full 
tests  are  underway. 

(4)  The  design  and  manufacture  of  the  equipment  for  the  experiment 
to  measure  the  effects  of  periodic  unsteadiness  on  the  aero- 
dynamic and  heat  transfer  performance  of  turbine  blades  has 
been  completed. 

(3)  Comparative  studies  between  Transonic  and  Subsonic  Turbines 
have  been  carried  out. 
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. AERODYNAMIC  AND  HEAT  TRANSFER  TESTING  OF  THE  FOUR  TRANSONIC  BLADE 
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PROFILES  AT  ZERO  INCIDENCE 

1 . 1 Blade  and  Cascade  Geometry 

The  convencional  transonic  blade  profile  Is  usually  designed  with 
a straight  suction  surface  with  the  throat  at  the  trailing  edge  [11],  as 
shown  in  Figure  1(a).  This  was  chosen  as  the  reference  blade.  The  other 
three  blade  profiles  are  variations  of  this  profile  to  a large  extent. 

Thus,  the  blade  shown  in  Figure  1(b)  has  an  unguided  expansion  on  the 
suction  side.  The  expansion  fan,  which  is  developed,  helps  to  weaken  the 

trailing  edge  shock  system.  As  a result  of  the  added  expansion,  the 

trailing  edge  dimension  comes  out  smaller  than  that  of  the  reference 
blade;  and  so  it  became  necessary  to  design  the  blade  shown  in  Figure  1(c), 

which  has  the  same  profile  as  the  reference  blade,  but  with  a thin  trailing 

edge.  The  last  profile  Figure  1(d)  shows  a longer  blade  with  a convergent- 
divergent  noizle.  A summary  of  the  blade  and  cascade  geometry  is  given 
below. 

1.1.1  Reference  Convergent  Transonic  Blade  with  Straight  Suction  Back 
aj Thick  Trailing  Edge  (Blade  1) 

Blade  chord  a = 66.0  mm 

Blade  spacing  Q = 49.5,  i.e.,  = 0.75 

Stagger  angle  y = 51° 

Inlet  flow  angle  = 60° 

Exit  angle  B,,  = 25° 

Trailing-edge  thickness  tc  = 2.8  mm,  tela  - 0.0424 
Number  of  blades  in  the  cascade  = 6 

with  Blade  '3  suction  side  and  Blade  4 pressure  side 
instrumented  with  stat  ic-prc-s.sure  taps. 
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1.1.^  Wiiolly  Convergent  lilade  witii  Stra:y;ht_  Suc_tijon  Hack  ]^ut  I’liin 
(railing  (Blade  2) 

Blade  chord  a - 64.2  mm 

Blade  spacing  g = 46  .1:>  mm,  i.e.,  g/c  = (J.7j 
Stagger  angle  y = 51° 

Inlet  flow  angle  = 6U°  , exit  angle  62  * 25° 

Trailing-edge  tnickncss,  te  = 1.3  mru,  te/c  = 0.02025 
•Junber  of  blades  in  the  cascade  • 6 

with  blade  3 suction  side  and  Blade  4 pressure  side  instrumented 
witli  static -pressure  taps. 


1.1.3  Transonic  Turbine  Blade  with  Ung^uide^  Zxpan^sioji  oti  the  _?uctioii^ 
Side_  (J’lug  fiozzle)  (Blade  3) 

Blade  chord  c = 06. 0 mn. 

Blade  spacing  g = 49.5  nm,  i.e.,  aln  = 0.73 
Stagger  angle  y “ -'i* 

Inlet  flow  angle  = b0°  , exit  angle  = 23' 

Taai  ling -edge  tiiickncss,  te  = 1.4^  r(m,  tele  - 0.0221 
.iumoer  of  ulades  in  tlie  cascade  = 0 

witii  nlade  i suction  sitie  and  Blade  4 pressure  side  instrumented 


wit.i  static-pressure  taps. 
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1.1.4  Blade  With  ConverRent-Divergent  Channel 
Blade  chord  e = 76.0  mm;  Throat  0 “ 18.30  mm 
Blade  spacing  g = 49.6  mm,  i.e.,  g/a  = 0.64 
Stagger  angle  y • 46.5° 

Inlet  flow  angle  = 6u°  , Exit  angle  ^2  “ ^5° 

Trailing-edge  thickness,  te  = 1.14,  i.e.,  te/a  = 0.015 
;;umber  of  blades  in  the  cascade  = 6 

with  Blade  3 suction  side  and  Blade  4 pressure  side  instrumented 
with  static  pressure  taps. 

1.2  Blade  Instrumentation  for  VKI  tunnel: 

The  blade  velocity  distribution  was  measured  at  the  mid-span  by 
pressure  taps  on  the  pressure  side  and  suction  side  of  two  neighboring 
blades  such  that  the  instrumented  blade  surfaces  formed  the  blade  passage. 

1.3  Inlet  Flow  Field 

Tile  Mach  number  distribution  at  the  cascade  inlet  was  derived  from 
wall  static  pressure  in  a plane  ''0.1c"  ahead  of  the  leading-edge  plane. 

Inlet  Macli  number  variation  as  a function  of  exit  Macli  numiier  was  measured 
and  tne  inlet-angle  variation  in  tiie  transverse  direction  was  also  measured. 
;,'o  influence  on  tne  outlet  Mach  number  was  observed. 

1 . 4 i.'.iscade  Flow 

Tile  blade  velocity  is  expressed  by  tiie  local  .‘iacli  number  calculated 


from  tne  local  static  pressure  on  tlie  Plade  and  t.ie  total  pressure  upstream 


r - ^ 


of  the  cascade.  Measurenents  were  taken  on  tlie  suction  and  pressure  sides 
of  the  blade.  Tlie  Mach  number  distributions  are  plotted  in  Kip.ures  2 
throup'.i  5. 

1.5  blade  Performance  as  measured  in  the  VKI  Wind  Tunnfil: 

The  blade  performance  was  evaluated  from  the  downstream  wake  traverses 
behind  the  blades  at  an  axial  distance  of  ''0.3c?"  behind  the  trailing  edge 
in  all  cases. 

1.5.1  Lo^£S 

The  variation  of  blade  losses  as  a function  of  the  outlet  Mach  number 
for  all  the  blades  is  platted  in  Figures  6 througli  9. 

For  Blade  1,  the  convergent  reference  blade  with  straiglit  suction  back, 
for  exit  Mach  number  up  to  =0.7  blade  losses,  ran  up  to  8%.  This  higii 
level  of  loss  was  apparently  caused  by  trailing  edge  thickness  (t<?/c?  = 4:.' 
as  compared  to  tele  = 2.025%  in  the  second  blade,  and  tele  = 2.2%  in  tae 
third  blade) . 

Local  supersonic  zones,  lambda  ("X'')  sliocks  and  diffusion  losses 

cause  a sudden  rise  in  the  level  of  tae  losses  from  M^  > 0.7  witli  maximum 

loss  of  about  11%  occurring  at  about  “ 0.85. 

The  decrease  in  the  level  of  losses  from  M.,  = 0.85  to  around  M,  

1.3  is  linked  with  the  fact  that  tlie  sliocks  became  more  oblitiue  and  flow 

reattachnent  occurred.  Tlie  new  increase  of  for  K,  > .1,  . up  to  tue 

2 2 ic£"-vt 

limit  loading  1.59,  is  due  to  the  increasing  strength  of  tae 
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left-running  trailing  edge  shock  which  caused  the  boundary  layer  to  separate 

without  reattachment  on  the  suction  side. 

Losses  recorded  for  M„  > M„  , . . . depend  on  the  total  pressure  losses 

2 2 Zt-mt 

due  to  the  blade  boundary  layer  ^ . the  trailing-edge  shocks 

and  the  mixing  process  between  the  trailing  edge  plane  and  the  measuring 

plane  as  well  as  the  outlet  velocity.  The  ^ remains  roughly  constant, 

while  (AP^)  , and  mixing  losses  increase  as  the  exit  Mach  number  increases. 

0 snook 

For  Blade  2,  (Figure  7),  the  convergent  blade  with  straight  suction 
back  and  thin  trailing-edge  thickness  {tela  = 2,025%)  for  exit  Mach  number 
up  to  “ 0.75,  losses  were  decreasing  until  they  reached  a value  of  4.7%. 
This  again  is  probably  due  to  flow  acceleration  and  thinning  of  the  boundary 
layer . 


As  in  Blade  1,  local  supersonic  zones,  lambda  ("A")  shocks,  and 

diffusion  losses  cause  a sudden  rise  in  the  level  of  the  losses  for 

M2  > 0.75  with  maximum  loss  of  about  6^  occurring  at  about  M2  ” 0.95. 

The  losses  then  decrease  as  in  tne  first  blade  from  = 0.95  to  around 

M . » 1.3  as  the  shock  becomes  oblique  and  possible  flow  reattachment 

2 aest-cm 

occurs  after  separation.  The  pattern  in  whicli  the  losses  decrease  from 

It  , . is  very  similar  to  that  of  Blade  1 (straight  suction  back  with 

2 aesvtjn 

thick  trailing  edge)  . 

Tlie  new  increase  in  u»  for  M-  > , . up  to  t'ne  limit-loading  • . .. 

2 2 aesim  ^ 2 ci'^nt 

is  due  both  to  mixing  losses  and  to  the  increasing  strength  of  tiie  left-running 
trailing  edge  shock,  which  causes  the  boundary  layer  to  separate  without 
reattaciiment  on  tne  suction  side.  The  very  sharp  increase  in  the  level  of 
the  losses  after  M.  ; • is  similar  to  that  of  blade  1;  at  around 


A 
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M2  • 1.5  each  reaching  a comparable  level  of  losses. 

For  blade  3 (Figure  8) , the  plug  nozzle  blade  with  expansion  on  the 
suction  side  (te/c  “ 2.21%),  for  exit  Mach  numbers  up  to  “ U.7,  losses 
were  decreasing  and  readied  their  minimum  value  of  5.5%,  probably  due  to 
flow  acceleration  and  thinning  of  the  boundary  layer. 

As  in  other  blades,  local  supersonic  zones,  lambda  ("X")  shocks  and 
flow  diffusion  cause  a sudden  rise  in  the  level  of  the  losses  for  M2  ^ 0.7, 
with  maximum  loss  of  about  7.5%  occurring  at  M2  “ O.V. 

The  losses  then  decrease  as  in  tne  first  and  second  blades  from 

M_  « 0.9  to  , . ■ 1.3  as  the  shock  becomes  oblique  and  possible  flow 

2 Zaesign 

reattacliroent  occurs  after  separation.  Tiie  decrease  is  fairly  sharp,  unlike 
Blades  1 and  2,  and  the  losses  flatten  out  between  M2  “1.2  and  M2  "1.4 
providing  a reasonable  operating  zone  witli  (u)  - 5%)  . This  region  of 
loss  lies  between  5%  and  5.3%. 

Losses  increase  for  M2  > limit-loading 

but  not  as  steep  as  in  the  first  and  second  blades.  This  is  probably  due 
to  the  comparatively  smaller  increase  of  strength  of  the  left-running 
trailing  edge  shock. 

For  Blade  4 (Figure  9),  ttie  convergent-divergent  blade,  {tc/c  = 1.5%), 
losses  were  first  decreasing  apparently  due  to  flow  acceleration  and  tninning 
of  the  boundary  layer.  Losses  reached  a low  value  of  b.5%  at  M2  “ 0.75 

As  the  pressure  ratio  increased,  local  sonic  bubbles,  lambda  C'X") 
sliocks,  and  flow  diffusion  increased  the  loss  level  to  a maximum  of  8.2j,» 
at  M , " O.80.  Losses  tlien  decreased,  as  in  Blades  1,  2 and  3,  but  at 
a very  sharp  rate  (from  M2  ” 0.86  to  around  M2  “ 

L 


1.3)  as  the 
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shock  became  oblique  (weaker),  and  as  flow  reaccaciimenc  occurred  after  separa- 
tion. A very  low  level  of  losses,  3.5%,  was  recorded  at  the  design  exit 
Mach  number  region  and  remained  low  (below  4%  between  - 1.2  to 
1.5),  providing  a good  working  range  of  low-loss  exit  Mach  number. 

A comparison  of  the  loss  curves  (Figure  4.19a)  shows  clearly  that 
each  profile  has  a superior  performance  in  a different  Mach  number  range 
and  presents  a unique  potential  when  cooling  problems  and  problems  of  structural 
integrity  are  carefully  examined  along  with  optimization  of  efficiency. 
Remembering  also  that  lower  blade  ctiord  a meatis  increased  secondary  losses 
tends  to  suggest  a more  broadly  based  assessment  of  losses  to  determine  an 
optimum  profile  for  a particular  application. 

1.5.2  Outlet  Angles 

The  outlet  angles  shown  on  the  performance  curves  were  measured  with 
reference  to  tiie  tangential  direction  (Figures  6 through  9) . 

For  blade  1 (Figure  6),  tlie  convergent  reference  blade  with  straigiit 
suction  back,  decreased  slightly  from  25.5“  to  24*  between  H2  = O.b  and 

S ' ^^2  design  ' ”2  ^ ^^Idesign’  increases 

slowly  up  to  - 1.4  in  a linear  fashion.  At  '’2  = loading^ 

the  deviation  was  about  10“ . 

For  Blade  2 (Hgure  7),  tiie  convergent  blade  witli  strai;,!it  suction  back 

but  tiiin  trailing  edge,  decreased  sliglitly  from  25.5“  to  23.5“  between 

= 0.5  and  .l^  = 1.1.  A gradual  increase  was  noticed  from  = 1.1 

to  . =1.3  wlien  increased  back  to  25“  . the  desij\n  exit  angle. 

2 lesvnn  2 
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For  Ml  ^ M,  tile  exit  angle  increased  very  sliarply  as  in  blade  1. 

For  Blade  3 (Figure  8),  the  plug  nozzle,  ,^2  decreased  from  28°  to 

8,  » 24.3°  between  “ 0.(>  and  !L  • 1.3.  For  M,  > M.,  . . , the  exit  angle 

z z z z z ^esion 

increased  very  sharply  as  in  blades  1 and  2 in  a linear  form.  At 

tL  “ M,..  “ 1.67  the  deviation  of  about  10°  was  again  recorded. 

2 llzmit  Load^nc 

For  Blade  4 (Figure  9),  the  exit  ar '.le  changes  very  little  for  all 

the  ranges  of  the  M up  to  M . . , fluctuating  by  about  0.3°  from 

22  aesvnn  ^ ’ 

a mean  of  24. j°.  (The  design  exit  angle  is  25°. ) After  M • (M  * 1.3), 

^ 2 ctesvcTK 

a sharp  increase  in  the  value  of  the  deviation  is  noticed,  already  amounting 
to  a 2*  deviation  at  .12  1.47,  which  suggests  a possible  large  deviation 

of  the  exit  angle  at  nigher  exit  Macii  numbers  similar  to  the  other  blades. 

As  siiown  in  Figure  10,  there  are  noticeable  differences  in  tiio  exit 
angle  behavior  for  each  blade  in  the  Mach  number  range  investigated.  The 
striking  difference  is  tiie  rather  large  deviation  in  blade  3,  close  to  3°, 
at  low  Mach  numbers.  This  is  not  surprising  since 

blade  3 is  the  plug  nozzle  witli  unguided  expansion.  At  these  low  subsonic 
Mach  numbers  it  acts  ns  an  inefficient  diffuser.  This  deviation  tiien 
disappears  as  the  design  supersonic  exit  Mach  number  is  approached. 

Anotlier  striking  point  is  that  a very  sharp  and  higli  deviation  is  noted 
in  all  the  cases  at  M,  > i.4,  i.e.,  just  a little  over  M- 
This  is  associated  witn  tue  over-expansion  at  the  trailing  cdRe  corner 
and  the  effect  is  even  seen  on  tue  wake  inclination. 
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1.5.3  Effect  of  Reynolds  Number  Variation 

In  addition  to  the  normal  test  series  performed  at  Reynolds  numbers 

between  10*  and  10*,  the  reference  blade  (Blade  1)  was  tested  at  constant 

inlet  total  pressure  of  P,,  “ P ^ 1000  inmllg,  constant  Reynolds  number 

01  atrmsp 

of  10^  (corresponding  to  tlie  highest  total  pressure  for  tne  hignest  Mach 
number  in  the  normal  test  series),  while  the  downstream  pressure  was  varied 
with  a back  pressure  valve. 

The  results  of  these  tests  showed  that  the  Reynolds  number  effect 
w->s  negligible  in  the  test  range  of  10*  to  10*.  No  noticeable  change 
was  seen  on  the  pressure  distribution.  Both  the  peak  position  and  its 
value  remained  practically  unchanged,  and  the  pattern  of  the  shock-system 
remained  exactly  the  same. 

Effect  of  Blade  Solldi.ty_qnJ^rf^rmance 
Tests  were  also  carried  out  to  determine  the  effect  of  stage  solidity 
on  cascade  performance.  Reference  Blade  1 was  tested  at  g/c  = 0.81 
and  y/t*  = 0.095  (former  tests  were  conducted  at  y/c?  = 0.75).  In  addition 
Blade  3 (plug)  was  tested  at  g/c  = 0.695. 

Careful  examination  of  the  blade  velocity  distribution  (compare 
Figure  3 aiid  Figure  11)  and  tae  scnlieren  photograpiis  show  a ciiange  in  tne 
locations  and  inclinations  of  tiie  left-running  shocks  for  the  same  values 
of  M.  . Also,  as  the  blade  spacing  increased,  tne  shock  got  weaker  and 

X’S 

some  flow  separations  resulting  in  tlie  sliock-boundary-layer  Interactions 


were  avoided. 
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1.5.5  Downstream  Wake 

Analysis  of  the  aowustreain  waite  using  tiie  downstream  traverse  taken 
at  0.3c?  behind  tne  blades  and  sciilieren  photos  shows  an  interesting  effect 
of  Mach  number  on  the  wake.  As  increases  the  flow  nonuniformity 
increases.  A characteristic  effect  is  the  increase  in  tiie  wake  depth 
during  the  transition  from  subsonic  to  supersonic  exit  Mach  number.  'Ihe 
width  of  the  wake  changes  also. 

For  the  different  blade  profiles,  tiie  effect  of  trailing  edge  thickness 
can  also  be  seen  clearly  in  the  relative  width  of  tlie  wake. 

Careful  examination  of  the  sciilieren  pliotos  of  tiie  wake  flow  shows 

that  it  consists,  under  certain  flow  conditions,  of  von  Karman  vortex 

streaks.  Because  of  tlie  long  exposure  time  of  tlie  camera,  an  estimate  ol 
the  sliedding  frequency  of  tiie  vortices  can  not  be  made  from  tlie  evaluation 
of  the  present  schlieren  pictures. 

Further  experiments  are  planned  in  conjunction  with  VkI  to  determine 
the  shedding  frequencies  oy  directly  measuring  the  pulsating,  wake  pressure 
with  a high-frequency-response  .iullte  pressure  transducer. 

1.6  Cascade  Performance  as  measured  In  the  MIT  Blowdown  K.uility: 

1.6.1  Blade  Pressure  Distribution  as  Measured  in  tlie  Mil  Blowdown 
I'ac  il  ity 

Plots  of  the  blade  pressure  (Macii  numi'er)  .1  i s t r i but  ions  are  sliowni  in 
Figures  12  and  13.  A very  fast  expansion  takes  p I ao<-  on  the  suction  sidi- 
witii  t iic  .‘lach  number  reac.iing.  its  first  peak,  at  a loi.it  ion  waere  tne  lirst 
surface  roucliness  (biimt'l  is  located  on  the  hi  n'e  . ihe  flow  leci-lerates 
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a bit  and  quickly  starts  to  accelerate  again,  until  it  liits  the  left-running 
shock,  from  the  neighboring  blade,  downstream  of  wiiich  tlie  flow  decelerates 
and  later  accelerates  towards  the  trailing  edge  on  the  suction  side. 

As  the  pressure  ratio  increases,  the  shock  moves  downwards  towards 
the  blade  trailing  edge  and  this  represents  a typical  pattern  of  pressure 
distribution  on  the  suction  side  for  all  the  blades. 

On  tlie  pressure  side,  a rapid  acceleration  occurs  downstream  of  the 
stagnation  point,  then  quickly  stops  close  to  the  leading  edge.  The  flow  then 
gradually  accelerates  to  its  maximum  value,  wliich  occurs  at  the  blade  trailing 
edge.  Once  tne  flow  is  ciioked  there  is  very  little  change  in  tlie  pressure  side 
pressure  distribution. 

All  tac  blade  profiles  tested  featured  the  same  typical  variation  in 
Mach  numoer  (pressure)  distribution,  differing  only  in  absolute  values  of 
tue  peaks  due  to  strength  and  locations  of  the  shocks  hitting  the  suction 
sides.  Since  all  the  blades  sliowed  tlie  same  tendency  to  peak  at  nearly  the 
same  place  on  tne  suction  surface,  it  would  seem  logical  to  attri- 
bute this  peak  to  tne  geometry  (large  curvature)  of  the  blade. 

1 . b . 2 Comparison  Between  Blade  Pressure  Distributions  Obtained  at  the 
M B lowdown  Cascade  Facility  and  tb.e  VKI  High  Speed  Wind  Tu nn e^l 

A detailed  ei'iiiparison  was  made  between  the  .Mach  number  distribution 
around  tlie  blade  obt.iined  in  the  MIT  hot  blowdown  cascade  facility  with 

I 

tu.it  ootaliied  in  the  Viil  iiigh  speed  windtunnel  (see  Figures  2 throng. i 5 
and  ligures  12  ami  13).  Jii  the  wiole,  tlie  overall  pattern  of  the  .Mach 
nuribi'r  distribution  w.is  similar.  l.isL  acceleration  from  tae  st.ignation  point 
l.u  suction  side  le  lu..  to  tai  tir.tt  peal  of  Mach  nuiimr.  hut  t "e 

location  anu  .lusnluti'  value  (nacnitude)  of  tne  pe.ik  was  different.  ^ 


I 
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The  peak  shifted  forward  in  the  tests  conducted  on  the  hot  blowdown 

cascade  facility.  There  arc  several  possible  reasons  for  tiiis.  First, 

it  must  be  recalled  that  the  experiment  in  the  windtunnel  at  VK'I  was  a 

cold  flow  investigation  (/  IT.  = 1)  with  a low  level  of  turoulence  0.7%. 

wiiereas  experiments  on  the  Ml'l  cascade  facility  were  hot  flow  experiments 

(inlet  stagnation  temperature  of  450®K,  I’  !'£,  = 1.5),  with  a high 

a r 

turbulence  level  (10%).  Scnlieren  photos  taken  at  VKI  showed  a distinct 
region  with  a separation  bubble  on  the  olade  suction  side.  A combined 
effect  of  large  blade  curvature  and  surface  roughness  (bump)  could  have  triggered 
this  separation  while  tiie  distinct  region  served  to  generate  the  necessary 
turbulence  in  the  shear  layer  to  reattach  the  flow.  Also,  the  higli  freestreani 
turbulence  could  allow  a quick  turbulent  flow  reattachment  and  thus  create 
an  earlier  transition. 

Another  possible  reason  could  be  that  there  are  fewer  pressure  taps 
on  the  blade  tested  at  the  MIT  hot  blowdown  cascade  facility.  Tliere  were 

twenty  five  pressure  taps  and  twenty  five  heat  transfer  gauges  evenly  dis- 
tributed all  around  a single  instrumented  blade  as  compared  to  instrumenting 
two  blades  (that  formed  tne  central  channel)  with  thirty  six  pressure  taps 
on  the  suction  side  and  twenty  two  on  the  pressure  side  in  the  tests  in 
tiie  VnI  windtunnel.  Thus  at  HIT,  pressure  taps  were  not  as  closely  located 
and  the  exact  location  of  the  peak  could  have  been  missed. 

Hut  on  the  whole  the  overall  pattern  of  the  pressure  distribution 
was  the  same.  Ttie  correct  peaks  of  tne  Macii  number  before,  at  and  alter 
the  shock  could  not  oc  exactly  located  but  t b.e  distribution  still  sliows 
the  existence  and  mt'Viment  of  the  sliocks  correctly. 

A good  agreement  is  seen  on  tin.  pressure  sine  witi.  graviual  nunbi  r 

increase  tip  to  tne  trailing,  edge. 

1 
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1.6.3  Blade  Heat  Transfer  Distribution 

Tlie  blade  beat  transfer  distributions  are  plotted  in  Figures  14  through 

19.  I'lots  of  the  local  Nusselt  number  .Nu  as  a function  of  the  relative 

— 0 8“ 
coordinate  X * X/a\  and  plots  of  Nu/Ke  ‘ versus  X are  all  siiown. 

(Additional  curves  of  .Nu/Re*^'^  versus  X and  hu/Re^'^*^  versus  X were  drawn 

for  the  blades.  (Figure  20.) 

To  aid  the  understanding  of  the  curves  for  iieat  transfer  and  static 
pressure  distributions,  a typical  pattern  of  heat  transfer  and  of  dimen- 
sionless pressure  coefficient  have  been  drawn  on  polar  diagrams  around  tiie 
blade  surface  (Figures  and  22). 

On  tlie  whole  the  external  flow  pressure  distribution  and  level  of 
turbulence  have,  as  expected,  greatly  Influenced  the  boundary  layer  flow 

which  in  turn  practically  dictates  the  heat  transfer  structure  on  the  blade 
Surface.  It  is  remarkable  to  see  a similar  effect  of  shock  boundary  layer 
interaction  on  both  the  pressure  and  heat  transfer  distributions  around 
the  blade,  may  be  seen  in  these  figures  also. 

Generally  the  pattern  is  as  follows.  On  the  suction  side  from  the 
leading  edge,  the  pressure  distribution  indicates  a very  fast  acceleration 
to  a minimum  pressure,  a condition  very  favorable  for  tiie  formation  of  a 
laminar  boundary  layer,  and  thus  the  <usselt  number  decreases.  Transition 
from  laminar  to  turbulent  flow  then  caused  a jump  in  the  Nusselt  number. 

The  Nusselt  number  tl ereaf ter  fluctuates  up  and  down  similar  to  the  beliavior 
of  the  pressure  distribution  in  this  area,  until  it  hits  the  left-running 
shock  coming  from  the  trailing  ec'ge  of  the  noigiiboring  tilaJe.  Downstream 
of  the  shock,  tlie  Nusselt  number  drops  significantly  and  later  starts  to 
increase  as  tne  flow  accelerates  toward  tne  trailing  edge. 

On  the  pressure  side,  from  the  li'ading  edge  the  pri'ssure  ilist  ribut  ion 
indicates  again  a rapid  acceleration  conducive  to  the  maintenance  of  .i  laminar 
boundary  layer,  and  .is  a result,  tue  Nusselt  number  falls.  Transition  oicurs 
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at  s/c  = 0.3,  and  tiie  iJusselt  number  increases  all  the  way  towards  the  trailing 
edge  due  to  the  rapid  decrease  in  pressure  up  to  the  trailing  edge  on  the 
pressure  side.  This  increase  in  the  Musselt  number  is  caused  by  the  thinning 
of  the  boundary  layer  due  to  favorable  pressure  gradient. 

It  must  be  mentioned  that  a hign  level  of  heat  transfer  was  recorded 
close  to  the  trailing  edge,  amounting  to  about  75%  of  the  mean  value 
measured  at  the  leading-edge  region.  It  was  also  noticed  that  zones  witn 

variation  in  heat  transfer  rate  generally  coincide  with  zones  with 
variation  in  pressure  distribution  (Mach  number). 

In  all  cases,  straight  lines  nave  been  used  to  connect  all  the  data 
points,  because  of  the  relative  distances  between  the  data  points,  tne 
true  picture  of  the  heat  transfer  distribution,  in  particular  the  exact 
values  and  locations  of  the  peak,  may  not  be  exactly  the  same  as  shown. 

The  heat  transfer  gauge  (3/32”  diameter)  itself  is  averaging,  over  a certain 

area,  across  which  the  .'iusselt  number  could  !iavc  varied  very  significantly. 

Thus  tlie  so-called  local  heat  transfer  distribution  is  in  fact 

an  average  over  a certain  small  area  in  that  region. 

1.6.4  Comparison  Between  Experimental  Heat  Transfer  Data  and  an 
Available  NASA  Theoretical  Prediction 

1 .6.4.1  Theoretical  Predict  ion  Based  on  Hot  Blowdown  Cascade 
blade  Pressure  Distribution  (as  Input) 

A comparison  between  the  oxper imi'ut a 1 results  and  values 
obtained  by  available  tlieoret  ical  method  is  shown  in  I'igures  2! 
througii  28.  A modified  NASA  Compressible  l.aminar  and  Turbulent 
>^1)11 'darv  Laver  Program  fl|  is  ussi  i"  iiredi.t  tin  Liai  transt.r 


d istr  i but  i on . 
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As  input  into  the  program,  the  experimental  pressure  distribution 
obtained  from  the  hot  blowdovm  cascade  tests  was  used.  The  curves 
(Figures  23  through  28)  did  not  agree  point  by  point  but  on  the  whole 
the  general  pattern  was  the  same,  and  it  was  an  agreement  good  enough 
for  such  a complex  quantity  as  the  heat  transfer  coefficient.  There 
are  obvious  reasons  for  the  differences  in  the  curves.  Surface  curva- 
ture, surface  roughness,  initial  turbulence  level  of  the  freestream 
and  shock  boundary  layer  interactions  are  not  taken  into  account  in 
NASA  program  except  insofar  as  they  affect  the  pressure  distribution. 

On  the  experimental  side,  the  physical  size  of  the  gauge  (3/32" 
diameter  x 0.020"  thick  aluminum  discs)  necessitates  an  average  measurement 
over  a certain  surface  area,  across  wliich  the  heat  transfer  could  have  varied. 

The  most  obvious  discrepancy  between  the  two  results  occurred  as 
a result  of  two  different  locations  for  the  transition  point.  The  experi- 
mental data  showed  an  earlier  transition  to  turbulent  flow  than  the 
theoretical  prediction. 

There  are  several  possible  reasons  for  this.  The  theoretical  prediction, 
as  mentioned  earlier,  does  not  take  into  account  tiie  freestream  turbulence, 
and  also  assumes  an  adiabatic  wall.  liie  experiment  on  tne  other  uand  aas 
a nigli  level  of  freestream  turbulence  wnich  tends  to  destabilize  the  laminar 
boundary  layer  and  forces  an  early  transition,  while  the  heat  transfer 
across  the  boundary  layer  to  the  blade  acts  at  the  sami'  time  to  staoiiize 
tile  iioundary  layer  by  uissipatinp  tne  energ-.’  of  turimlence . Tiiis  will  tend 
to  delay  transition.  The  net  effect  wlien  coupled  witu  Llie  effect  of  tne 
lilade  surface  rmigimess  (the  blades  in  fact  di'  have  bumps  in  tiiis  region 
of  interest)  could  iiave  caused  tiie  earlv  transition  . 
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On  the  pressure  side  however,  the  experimental  data  on  the  heat 
transfer  suggests  that  instead  of  having  a transition  point,  a transition 
region  could  in  fact  have  existed.  This  is  possible  under  the  favorable 
pressure  gradient  in  this  region. 

1 . 6 . A . 2 Theoret ical  Prediction  Using  VKI  Wind  Tunnel  Cold  Flow 
blade  Pressure  Distribution 

Another  series  of  theoretical  predictions  were  done,  this  time 
using  the  experimental  pressure  distribution  obtained  from  cold  flow 
VKI  wind  tunnel  tests. 

The  three  curves  of  heat  transfer  distribution  (experimental 
data,  theoretical  prediction  based  on  the  high  turbulence  hot  flow 
blade  pressure  distribution,  and  theoretical  prediction  using  the  low 
turbulence  cold  flow  pressure  distribution)  are  shown  in  Figures 
2^  through  31. 

The  prediction  based  on  low  turbulence  cold  flow  is  very  similar 
to  that  based  on  high-turbulence,  hot-flow  except  for  the  location  of 
the  suction  side  transition  point.  There  are  also  differences  in  the 
absolute  value  and  location  of  peaks  but  apart  form  this,  there  is  very 
good  agreement.  The  low-turbulence,  cold-flow  has  the  most  delayed 
transition  as  expected.  This  is  definitely  the  effect  of  tne  difference 
in  the  freestream  turbulence  on  the  transition  of  the  boundary  layer. 

1.6.5  Transition  Point 

Sever.il  investigations  [2,  3,  * and  j ] n.ave  sliown  tnat  tne  t;ain 
flow  level  of  tnrimlence  exerts  considcrahli-  influence  on  the  stability 
of  a l.vninar  boundary  layer  ami  on  aeat  transfer.  LI  f«  i t of  turbtm.ice  (>n 
lo.i.ii  s in  tur'.iir.e  c.ascades  is  ..lniwn  It-  tiv  stiuct  ire  of  t.ii  tnnn.i.arv  l.i.er. 
tri  ti  in  lossc-  am'  int.-nsitv  of  .eat  exch.an  e an  ..Iter  •...tl.d'v,  ■ .ers.  . 


-19- 


depending  on  whether  the  flow  regime  in  the  boundary  layer  is  laminar  or 
turbulent.  It  is  therefore  essential  to  know  the  true  location  of  transition. 

normally,  in  the  absence  of  high  scale  surface  roughness,  flow  in  the 
boundary  layer  changes  from  laminar  to  turbulent  either  through  instability 
resulting  from  the  growth  in  tnickness  of  tlie  laminar  boundary  layer  or 
through  turbulent  reattaclunent  after  laminar  separation.  Studies  on  the 
behavior  of  turbulent  boundary  layers  in  accelerating  flows  [6,  7]  have 
shown  that  when  the  acceleration  is  severe  enough  and  acceleration  parameter 
3 

A “ — > 3 X 10  , the  originally  turbulent  layer  undergoes  a reversion 

uX 

e 

towards  laminar  just  as  in  this  experiment. 

Unfortunately,  quantitative  data  on  the  influence  of  freestream  turbulence 
level  on  transition,  lieat  transfer  [2,  3]  and  losses  arc  scanty  and  tne  few 
references  that  are  available  substantially  differ. 

Zysina-llolozhen  [4]  and  others  have  been  studying  tlic  effect  of  tur- 
bulence on  transition  in  the  boundary  layer  of  gas  turbine  blades.  They 
came  up  witn  an  empirical  formula  for  calculating  tne  Reynolds  number  at 
the  beginning  of  transition  as  a function  of  the  level  of  turbulence 

Z , wall  to  gas  temperature  ratio  and  Mach  number: 


H = 0.71  10 

extr 


where  is  the  llach  number  of  the  flow  at  the  naxinun  velocity  region  before 

transition;  !■!.  is  the  Mach  numiier  of  tiie  flow  at  cascade  inlet;  and  H = 

0 

Plotting  H [A]  versus  7 fturbulenre  level)  for  tvpical  values 


1 


of  "f-,  rK.  and  ’vi  has  shown  that  R , decreases  rapidly  when  T increases 
from  1.5/«  to  3™,  indicating  tnat  the  transition  point  moves  forward  and  tnen 
remains  at  a relatively  fixed  position. 

Zysina-Moloznen  and  Kunst  [A]  noticed,  however,  that  far  > 4.51., 
the  dependence  alters  and  some  stabilization  is  observed  in  whicn  naving 

reached  tiie  value  of  “ 1-2  x 10^,  noticeably  ceases  to  decline  witii 

a further  increase  in  turbulence.  It  is  as  tiiougii  it  reaches  the  maximum 
possible  displacement  of  the  transitional  point  upstream. 

but  the  quantity  /I’^^^^cannot  truly  be  a good  measure  of  the  transition 
point  since  it  takes  no  account  of  acceleration  of  tiie  flow.  Re^  would 
probably  have  been  a muca  better  measure. 

In  tlie  experiment,  transition  is  controlled  by  tlie  streamwise  pressure 
liistory  (turbulent  reattacament  after  laminar  separation)  . The  rapid  rise 
in  Nusselt  number  wliich  occurs  at  xta  = 0.2U  on  the  suction  surface  and 
xtc  ~ 0.35  on  the  pressure  surface  is  taken  as  evidence  of  transition. 

hdiowing  the  e::act  location  of  transition  and  forcing  the  theoretical 
program  to  start  to  calculate  turbulent  i’oiindary  layer  from  there  would 
probably  greatly  improve  the  quality  of  the  theoretical  lieat  transfer 
prediction . 

1 . 6 . h er_ the  Blade  Leading  Kd^e  Keg i on 

It  has  been  sliown  by  numerous  investigators  that  with  increasing  flow 
tiirinilence , neat  transfer  is  intensified  over  tiie  entire  leading  surfaces 
of  cylinders  or  spiieres,  particularly  in  tne  vicinity  of  tlie  frontal  st.ig- 


nation  point. 
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Kestin  [2]  and  others  have  correlated  heat  transfer  in  the  frontal 
stagnation  point  of  a cylinder  as 


This  relationship  is  based  on  the  assumption  that  lieat  transfer  is 
intensified  in  this  case  by  eddies  generated  in  the  boundary  layer  near  tiie 
stagnation  point  by  flow  turbulence,  with  the  distance  between  tlie  eddies 
(wavelength)  being  inversely  proportional  to  the  Reynolds  number. 

Uyban  e_t  al_^  [8]  suggest  from  general  consideration  of  similarity 
theory,  tue  use  of  the  turbulent  Reynolds  number  (T,  ‘Re)  and  correlate 

ti 

data  by  tiie  expression 


.'^Uji  * -'Ju^  * j (Re_  ) aJu™ 

u u u u u 

A comparison  is  nade  octween  tlie  experimental  neat 
at  tiie  leading  edge  with  tiredictiona  based  on  empirical 
formulation  of  other  works  listed  below. 


(1.2) 

transfer  data  obtained 
and  semi -empirical 


1.6. 6.1  Kestln  and  Wood  Correlation 
Kestin  and  i.'ood  (2]  formulated  tliat  in  tiie  range  0 < 
heat  transfer  at  the  leading  edge  could  be  expressed  as 


1/2 

(T  Re  ' ) < 40, 
u 


(T,  •Re^^'^\ 

!t 

= 0.943  + 1.48  U- 

- 3.9h  1 

' u 

\ luo 

\c  ' ^ ' 100  / 
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In  liquation  (3)  the  turbulence  level  is  expressed  as  an  absolute 

fraction  (as  opposed  to  a percentage)  , and  the  Reynolds  number  is  based  on 

leading-edge  diameter  and  inlet  flow  stagnation  parameters.  All  the  tests 

were  carried  out  between  lb . 3 < T < 27.6. 

— u 

1.6. 6. 2  Smith  and  Kuethe’s  Correlation 

On  tile  basis  of  a semi-empirical  theory.  Smith  and  Kuethe  [y]  were  the 
first  to  suggest  tlie  use  of  a single  correlation  parameter  T 

u 

.\ccording  to  tlieir  theory 

= U.3762  + 0.0133  - 1 . 32 • lo"^(2^ •Re^'^")  (1.4) 


1 . 6 . 6 . 3  Mujumdar  and  Douglas's  Correlation 

Mujumdar  and  Douglas  [10]  observed  that  a better  correlation  could  ue 

made  in  terms  of  turbulent  Reynolds  numi<er,  Re  .,(^T  *Re)  , instead  of  tne 

..  ;< 

1/2 

single  parameter  2’^ ‘Re  , and  gave  tiie  following  correlation 

iiu  -4  _n  , 

"Y/2^  * U.ydl  + 1.017*10  Ke„,  + 2.74x10  Re^,  (1.5) 

Re  ■' 

1 . 6 . 6 . 4  Uyban,  Ep i k and  Kozlova's  Correlation 

Judging  from  tne  results  of  tneir  wide  range  of  experimental  investi- 
gation.i  on  neat  transfer  in  tiie  vicinity  of  tin-  front  tagn.ition  iioint  of 
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a cylinder  in  traverse  flow,  Dyban  e_c  al . [b]  concluded  that  freestream 
turbulence  affects  heat  transfer  in  tl»at  critical  point  in  the  same 
manner  as  does  the  Reynolds  number. 

Lsinf.  turbulent  Reynolds  number  rj«Ke,  tliey  came  out  with  the 
empirical  formulation 


C 

o 


.MU, 


-0 

u 


0.8  r'*Re 

_ _K_  

iso'o  + f '■’•‘Re 


(1.6) 


where 


and 


Ru^  is  .lusselt  number  at  stagnation  point  in  turbulent  flow 

Nuj,  is  the  Mussclt  number  at  stagnation  point  in  zero  turbulence 
u 

e is  the  ratio  of  the  two  values  above 
o 

is  the  percentapc  level  of  tur!)ulence 
Re  is  tne  Reynolds  number  based  on  leading-edge  diameter  and  inlet 
flow  conditions. 


iVhen  experimental  data  are  compared  wita  tnese  f ornulations , as  shown 
in  Figure  32,  the  results  snow  that  the  measured  Nu  . at  tlio  blade  leading 
edge  were  consistently  lower  by  as  much  as  30%  tlian  the  predicted  values 
in  all  cases  except  for  tlie  Smitn-Kuethe  correlation  in  vliiru  tne  differences 
decrease  to  less  than  7... 

There  is,  nowev'er,  some  question  about  the  Smith-Kuetne  correlation 
for  ^ 20.  As  Kuitne  himself  has  indicated,  tlie  theory  would  require 

1/2 

modifications  for  > 20.  Since  tin-  experiments  reported  here  were 

1 /2 

in  tne  ranee  In.  3 < 1'  *"  27. l.,  tiu’  ai-reenent  in'tween  tne  I'xper  iriental 


data  ana  tne  Smi  tn-knet  ,u  correlation  siionld  bo  voiwed  with  re.serv.n  ion . 
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Tliere  is  no  doubt  tiiat  tiie  measured  Nu  was  consistently  lower  tnan 
its  true  value.  An  explanation  for  this  is  tiie  possibility  tiiat  the 
leading -edge  gauge  could  not  have  been  placed  exactly  at  the  stagnation 
point  and,  even  if  it  were,  it  was  averaging  heat  transfer  over  an  area 
in  the  leading  edge  zone,  where  tiiere  is  a very  sharp  gradient  of  pressure 
and  iieat  transfer.  The  ratio  of  gauge  diameter  to  leading  edge  diameter  “ U.3. 

It  is  therefore  not  unreasonable  to  indicate  that  tlie  maximum  iieat 
transfer  at  tne  leading  edge  may  not  nave  been  recorded  but,  instead, 
averaged  out  over  a certain  area  of  tne  leading  edge  region. 

1.6.7  Effect  of  Mach  Number  on  Blade  Nusselt  Number 

j\r,  attempt  was  made  to  correlate  the  effect  of  liacii  number  on  tiie  blade 

«l 

heat  transfer.  Variations  of  hu/Re",  for  various  values  of  n (n  ■ d.3,  d.oo, 

0.8  and  1.0)  were  investigated  and  some  are  shown  in  figures  17  tnrougii 
20.  Stanton  number  ,as  a function  of  the  dimensionless  ciiordwise  location 

was  also  investigated.  Unfortunately,  tlie  shock  movement  on  tne  suction 
side  makes  it  ini|)ossible  to  arrive  at  a meaningful  quantitative  correlation. 

A plot  of  mean  .Nu/Xe**'^,  obtained  from  tne  area  measurements  of  tne 
Iieat  transfer  dist r lout  ion  curves,  for  tne  entire  olade  surface  is  plotted 

as  a function  of  exit  M.icli  iiumoer  in  Figure  13.  The  plot  shows  tiiat  t.ie  , 

mean  .leat  transfer  to  tlie  blade  increases  witii  Macli  mimiior. 

A comparison  is  also  made  between  the  me.n.ured  olade  mean  ..usselt 
number  v,it.i  several  available  experimental  d.it.i.  As  s.iowii  in  Figure 
our  experimental  d.ita  lie  in  t;ic  upper  reeion  (iiign  valuis)  whicn  is 


eXfie..  tell  as  a result  of  the  liigii  level  ul  turaiileiui  at  w.jI,.. 


epi  1 a 


Finally  a comparison  of  the  mean  heat  transfer  (Stanton  number) 
as  a function  of  exit  Macn  number  is  made  of  all  the  blades  tested. 

This  comparison  (Figure  35)  shows  that  the  average  heat  transfer  coefficient 
of  Blades  1 and  3 is  about  the  same  whereas  the  convergent-divergent  blade 
has  a iiigher  Stanton  number.  This  high  value  of  the  Stanton  number  is 
probably  due  to  the  higher  level  of  pu  over  the  surface  of  the  convergent- 


divergent  blade. 


II.  UPDATIKG  THE  INSTRUMENTATION 


In  updating  the  instrumentation,  a look  was  taken  at  the  effect  of 
the  temperature  difference  between  the  copper  slug  and  the  rest  of  the 
surface  of  the  blade  which  results  in  the  pulsing  of  the  thermal  boundary 
layer.  In  order  to  estimate  this  effect,  an  approximate  integral  method 
was  used  as  shown  below: 


Figure  (a) 


I +Ax=l  +I.+dq 

X X 6 

i . e . 1 + Ax  - I = 1 - + dc; 

X X ' 
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Thus  Equation  (2.1)  becomes 

6. 


dx 


pV  c Tdy  - c T ~ 
X p p <»  dx 


pV  dy  = 

X 


- k(|^) 

y=0 


Assuming  constant  p and  c we 


1- 

f 1 

dx 

V^(T^-T)dy 

J 

0 

have 


,3T, 
9y  0 


(2.2) 


Assume  that  the  momentum  and  temperature  profiles  are  similar  and  use  the 
7th  power  profile  to  obtain 


V 

X 

V 


(y/6) 


1/7 


(2.3) 


T 

T 


00 


-T 


0 

0 


(y/<5.j) 


1/7 


(y/fi)^^^-(6/6^)^^'  = e 


(2.4) 


where  Tq  is  the  wall  temperature  and  6 and  6^  are  the  momentum  and  thermal 
boundary  layers  respective! ' 

Using  Equations  (2.3)  and  (2.4)  in  (2.2)  gives 


I(i,/6)*"6V„(T,-T„)l 


(il) 


= (T  -T,J(6/^.J 


.■  nJ/7  3 


0 


‘0''  T'  3v  'V 


C~' 


1/7  '0 

(7/6.,)"^ 

V I IJ 


(2.5) 


(2.6) 
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Equation  (2.6)  into  (2.5)  gives 


1 d ,,6/7^.  ‘O  1 1 


(2.7) 


where  ^ - 6^/6  and  Tq  is  the  friction  coefficient  at  wall, 
For  turbulent  flow; 


V X -1/5 


6/x  = 0.37( ) 

V 


(2.8) 


T_  V 6 -1/A 

— ^ = 0.0225(-^) 
pV^ 


(2.9) 


Using  Equations  (2.8)  and  (2.9)  in  (2.7)  gives 


dx^^  ^ Pr  ^ 

Equation  (2.10)  may  be  rearranged  as 

,8/7  A/5  1/8  d ,^8/7  A/5,  5.62  -1/10 

(C  X ) ^ ) = -JT  ^ 


(2.10) 


(2.11) 


Solving  Equation  (2.11)  gives 


q = 


6t 

T 


= A . 56  Pr 


-7/9 


(x^,/x) 


9/10 


7/9 


(2.12) 


where  Xq  is  an  unstarted  heating  length  whicfi  represents  the  distance  from 
the  leading  edge  to  tlie  gauge  in  que.-tlon.  f((iinti<  (2.12)  in  (2.h)  gives 


-29- 


= ^-0233  Rex°-®  - (x^/x) 


(2.13) 


Equatior  (2.13)  gives  the  correct  power  for  the  Reynolds  number  for 
Xq  = 0 but  gives  a lower  dependence  on  Prandtl  number. 


Figure  (b) 


Thus  if  it  is  required  to  find  the  effect  of  the  first  gauge  on  the  second 
gauge  then  Xp  = 0.25";  x = 0.5";  x = 0.35".  Hence, 

1 Ur. 


= 


AT 


T„-1V 


(R1-R2: 


where  Aq  is  the  error  in  recorded  he.it  transfer,  and  R = [1  - (x^/x) ° ^ ^ ^ ^ ^ 
and  AT  is  the  temperature  difference  between  the  metal  slug  and  the  blade 
surface  (Ficc  . 36).  For  the  conditions  of  tb.i'  experiment  at  t = 1 sec., 

= 47‘’F  and  = 270“i.  Then  5q/d  = 0.0113. 

The  effect  of  tlie  fust  gauge  on  the  third  g.uige  is  even  smaller  (0.151). 
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Thus  it  was  reconunended  to  stagger  the  heat  transfer  gauges  as  shown  in 
Figure  (b)  in  order  to  keep  these  effects  to  less  than  1%.  The  obvious 
weaknesses  of  the  analysis  include  the  assumption  of  constant  density  and 
the  fact  that  the  effects  of  curvature  have  not  been  taken  into  account. 
However,  staggering  the  gauges  is  expected  to  keep  the  errors  very  low. 

The  photograph  of  the  new  meliumeuted  blades  are  shown  in  figure 
37.  The  static  pressure  tubes  are  embedded  further  inside  the  blade  so 
that  less  distortion  of  the  flow  over  the  blade  surface  is  produced. 

The  bridges  for  matching  the  heat  transfer  gauges  and  the  pressure  taps 
with  the  oscillograph  have  also  been  made  a permanent  feature  of  the 
facility.  The  master  for  the  manufacture  of  the  blades  was  smoothened  so 
that  the  lumps  on  the  blade  were  eliminated.  The  resultant  blade  profile 
had  a slightly  slower  expansion. 


1 

1 

i 

' 

1 
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III.  OFF-DESIGN  PERFORMANCE  OF  THE  REFERENCE  TRANSONIC  BLADE  PROFILE; 


The  off-design  performance  of  transonic  blades  is  of  great  importance 
to  the  designer.  Thus,  in  this  section  of  the  study,  an  investigation  of  the 
aerodynamic  and  heat  transfer  performance  of  the  reference  blade  profile 
figure  1(a)  was  carried  out  at  angles  of  incidence  other  than  0.  llie  tests 
have  not  yet  been  completed,  but  some  preliminary  results  at  one  exit  Mach 
number  of  1.14  has  been  obtained.  Thus,  while  no  final  conclusion  will 
be  drawn,  some  preliminary  comments  will  be  made. 

Figure  38  (a) , (b) , and  (c)  show  the  isotropic  Mach  number  distribution 
obtained  for  angles  of  incidence  equal  to  zero,  +5°  and  +10°  respectively 
for  an  exit  Mach  number  of  1.14.  The  movement  of  the  stagnation  point 
is  apparent,  and  it  appears  as  if  the  transition  point  extends  into  a 
transition  zone  which  seems  to  get  larger  with  increasing  angle  of  incidence. 
The  Mach  number  at  transition  also  appears  to  be  higher. 

Figures  39(a),  (b) , and  (c)  show  the  Nusselt  number  distribution  for 
angles  of  incidence  equal  to  zero,  +5°  and  +10°  respectively  for  an  exit 
Mach  number  of  1.14.  A comparison  of  the  distributions  shows  the  movement 
of  the  stagnation  point.  These  results  are  preliminary  and  further  testing 
which  is  underway  is  required  before  any  trends  can  be  recognized  and 


conclusions  drawn. 


IV.  EFFECT  OF  UNSTEADINESS  ON  THE  AERODYNAMIC  AND  HEAT  TRANSFER 


? 


terformance  of  turbine  blades 

Flows  in  gas  turbines  are  necessarily  unsteady.  The  unsteadiness 
arises  from  factors  such  as  blade  passage,  the  propagation  of  azimuthal 
nonuniformities  and  the  propagation  of  periodic  disturbances  which  originate 
upstream  of  the  turbine. 

These  disturbances  are  likely  to  have  a marked  effect  on  the  flow 
field,  the  base  flow  and  the  trailing  edge  shock  system  behavior  as  well 
as  on  the  heat  transfer.  The  aerodynamic  flow  around  the  trailing  edge 
of  the  blade  is  of  great  Interest  because  of  its  effect  on  blade  row 
losses . 

When  a transonic  cascade  is  operated  at  conditions  other  than 
limit  loading,  the  left  running  trailing  edge  shock  of  the  first  blade 
usually  interacts  with  the  boundary  layer  of  the  suction  side  of  the 
adjacent  blade  and  this  interaction  could  result  in  local  or  full  flow 
separation,  which  has  a great  influence  on  aerodynamic  losses.  Since 
the  rotor  typically  operates  at  a high  angular  velocity  with  many 
individual  blades,  a stator  could  see  in  tiie  neighbortiood  of  10^ 
fluctuations  per  minute.  Tiie  purpi>se  I'f  tliis  experimental  investigation 
is  to  examine  the  effects  of  tiie  high  frequency  periodic  unsteadiness 
on  the  performance  of  tiie  given  turbine  blade.  The  experiment  is  being 
designed  to  investigate  tlie  loliowing  pheni'mon.i: 

(i)  tile  aerodynami  and  he.it  Iroislio  [ roperties  in  t lie 
neighborhood  ot  the  I'sc  i 1 1 .1 1 i ng  shock-bound.iry  l.iyer 
interaction  im  tin  sue  t ion  -.ide  ot  the  adj.teent  blade; 


( i i ) t tie  I 1 f I . t " t 'he  I'l  1 1 1 1 o 


lie , 


• ,1.1 i 1 s t II  I'  111.  e 'll  t be 
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potential  flow  in  the  blade  passage; 

(iii)  the  effect  of  the  disturbances  on  the  base  pressure  behind 
the  stationary  blade; 

(iv)  the  accumulated  losses  in  the  wake. 

The  above  experiment  is  to  be  conducted  utilizing  the  MIT  linear 
cascade  tunnel  which  is  attached  to  the  hot  blowdown  facility  in  the  Gas 
Turbine  Laboratory.  In  the  experimental  investigation,  a disturbance 
generated  by  rotating  a small  elliptical  body  downstream  of  the  cascade 
blade  will  be  used  to  simulate  periodic  unsteady  effects.  Blade,  cascade 
wall,  and  traversing  rake  instrumentation  will  be  used  in  measuring  the 
aerodynamic  and  heat  transfer  properties  for  different  operating  conditions. 
A traversing  probe  will  be  used  downstream  of  the  cascade  to  measure 
cascade  losses. 

The  experiment  sliould  help  provide  a basic  understanding  of  the 
losses  which  are  inherent  in  all  gas  turbine  generators  which  utilize 
transonic  rotating  machinery.  It  should  also  better  help  to  understand 
the  cooling  requirements  for  this  type  of  turbine  blading. 

A small  air  turbine  was  chosen  to  provide  the  required  flexibility 
of  rotating  the  elliptical  piece  downstream  of  the  turbine  cascade.  The 
air  turbine  has  a variable  speed  of  from  0 to  lOO.OOORPM  while  operating 
off  the  shop  air.  The  company  which  was  riven  the  contract  for  the  air 
turbine  is  the  Barbour  Stockwell  Comnanv  of  Cambridee.  MasRarbuBett*! 

"’he  air  turbine  was  supposed  to  be  delivered  in  May  but  during  the  acceptance 
tests  problems  occured  with  the  bearings  of  the  turbine  abc've  90,000  RI’M 
and  delivery  was  not  possible  in  Mav.  It  was  later  decided  t('  run  the 
turbine  below  90,000  RPM  and  so  the  air  turbine  will  be  delivered  within 


two  weeks  of  the  time  of  writing.  Meanwhile,  components  of  the  air  tur- 
bine control  system  consisting  of  a magnetic  pickup  attached  to  the  air 
turbine,  an  electrically  operated  throttle  valve  and  the  air  cleaning 
system  have  been  installed. 

Work  is  also  underway  in  the  area  of  an  analytical  study  of  the 
effect  of  unsteadiness  on  the  aerodynamic  and  heat  transfer  performance 


of  the  blades. 
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V.  COMPARATIVE  STUDIES  BETWEEN  TRANSONIC  AND  SUBSONIC  TURBINES: 

For  this  comparative  study,  let  us  consider  a case  where  one  highly 
loaded  transonic  turbine  stage  is  used  to  replace  two  conventional  subsonic 
turbine  stages. 

It  is  assumed  that  the  mass  flow  rate  mg,  turbine  inlet  stagnation 
temperature  T^^^  and  pressure  Pqj^  > work-output  and  tip  speed  are  the  same 
in  both  cases.  In  addition,  the  axial  velocity  is  assumed  constant  in 
both  cases. 

For  the  Transonic  Turbine 
Number  of  stages  = 1 
Stage  Stagnation  Pressure  Ratio  = 4 
Turbine  Inlet  Stagnation  Temperature  = 1800°K 
For  the  Subsonic  Turbine: - 
Number  of  stages  = 2 
1st  Stage:-  Stage  Pressure  Ratio  = 2 

Turbine  Inlet  Stagnation  Temperature  = 1800‘’K 
2nd  Stage:-  Stagnation  Pressure  Ratio  = 2 


5.1  Casdynamic  Analysis  of  the  ^wo  Turbines 

Gasdynamic  analysis  of  the  Transonic  Turbine  Stage  has  been  carried 
out  and  the  results  given  in  Table  Number  1.  A similar  analysis  is  now 
carried  out  here  for  the  two  stages  of  the  Subsonic  Turbines. 

Turbine  l'*’^  Stage:- 
^01  " 


0 i 


0 1 


.-I 


.'■'S 

P 

I)  ) 


1 . 1 V) 
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•••  "03  = 


AT  = 1800  - 1560°K  = 2A0°K 
o 


specific  Work  Output 


AW  = cp  . ATq  = S ^^0 

2 2 

= 299154.93  m /sec 


For  the  2nd  stage. 


Toi  = 


01 

D 

03 


2, 


01 

''03 


= 1.155 


Toi  = 


ATy  = 210°K 


The  detailed  multi-stage  analysis  is  shown  in  Table  11. 


-37- 


Table  II 


Stage  Parameter 


1st  Stage 


2nd  Stage 


1 . Design  Parameters 


(a) 

Turbine  Inlet  Total 
Temperature 

1800°K 

1560°K 

(b) 

Stage  Stagnation 
Pressure  Ratio 

2 

2 

(c) 

Blade  Speed  (C) 

550  m/s 

550  m.s 

(d) 

Target  Total  to  total 
polytropic  efficiency 

0.9 

0.9 

2. 

Stagnation  Temperature  Drop 
Across  Stage 

= T - T 

0 01  03 

240°K 

210°K 

3. 

Stage  Specific  Work  Output 

AW  = C AT„ 
P 0 

299154,9 

260887. 

4. 


5. 


h. 


Change  In  Tangential  Vel 

Aw 

^Cy  = — 543.92  474.34 


Tangential  Components  of 
the  Velocities 


W = i^£iL±±) 

y'i  ^ 2 ' 

546.47 

512.17 

C = W . - U 
y3  y3 

-3.53 

-37.83 

W - = C - 
y2  y3 

-3.53 

-37.83 

C „ = W - + U 
y2  y2 

-3.53 

-37.83 

Rotor  Inlet  and  Kxit 

Velocities 


/c^  + c" 


v2 


657 .16  ro/s 


6 2 H ni  / s 


Tab  If  11  (tumt  inued  ) 


8. 


Stage  Parameter 


1st  Stage 


2nd  Stage 


'2  2 

W_  = /W  + W „ 

1 r X V2 

365.02  m/s 

367  m/s 

s=  A'  "A 

3 / w y3 

365.02  m/s 

367  m/s 

' 2 2 

W-  = /W  + W - 
3 / X y3 

657 . 16  m/s 

628  m/s 

Static  Temperatures  & 

Acoustic  Velocities  at 

Rotor  Inlet  and  Exit 

(a) 

c2 

1626. 19°K 

1401  .27°) 

(b) 

a =/y  - T 
2 / * m 2 

778.51 

722.67 

(c) 

c2 

T = T - 3 

1506.37 

1295.79 

(d) 

33  =/y  ^ T 
J / m J 

749.28 

694.94 

Mach  Number  at  Rotor 

Inlet  & Exit 

(a) 

Absolute  Inlet  Mach  No. 

M _ = C./A„ 
c2  2 2 

0.84 

0.87 

(b) 

Relative  Inlet  Mach  No. 

\2  = 

0.47 

0.51 

(c) 

Axial  Inlet  Mach  No. 

M , = C /a_ 
cx2  X 2 

0.47 

0.  Dl 

(d) 

Absolute  Exit  Mach  Ni' . 

M = c /a 

c3  3 3 

0.49 

0.53 

(e) 

Relative  Exit  Mai-li  No. 

M = W /a 

w3  3 3 

>'.8H 

;.')0 

J 


1 


Table  11  ( 
Stage  Parameter 

3 3- 

icontinued) 

1st  Stage 

2nd  Stage 

(f) 

Axial  Exit  Mach  No. 

M n = c /a 
cx3  X 3 

0.49 

0.53 

9. 

Stagnation  Temperature  at 
Rotor  Inlet 

Absolute  ^02^^ 

1800“K 

1560“K 

rp  1 

Relative  Tq2 

1679.81 

1455. 48°K 

10. 

Stagnation  Temperature  at 
Rotor  Exit 

Absolute  ^03^^ 

1560°K 

1350°K 

Relative  ^03^^^^ 

1679.81 

1455. 48°K 

11. 

Stage  Loading  Coefficient 

, C AT 

= p o 

0.99 

0.86 

12. 

2 

y 

Stage  Flow  Coefficient 
<p  = ^x 

0.66 

0.66 

y 

5.2  EsLlmato  of  the  Coolant  Flow  Regulrenieiits  Film  Cooling 

As  mentioned  above  in  Section  5.1,  film  cooling  of  varying  degrees 

oi  sopbisitication  offers  a more  practical  solution  for  high-temperature 

* 

turbine  blade  cooling. 


Film-cooling  effectiveness  data  on  flat  plates  gathered  at  M.I.T. 
Gas  Turbine  Laboratory,  which  Incorporates  the  effect  of  geometry  on 
film-cooling  effectiveness,  shows  that  for  double  row  of  holes,  with  in- 
jection angle  of  20°,  gives 


n.  , = 156.17 

isoth 

.24 

-1.35  „-0.68  X 

(m  0 ) (p  . 

-.24 

(5.1) 

= 156.17 

1 

1 

.24 

V - . 24 

-1.35  X M 

(Ky  ) (p  . y 1 

U 

u pc 

where  m = -7;—  ; 


7 


T 

c exit 
T 

q 


pu“ 


pU- 

R 


K 


0 


m/O 


M 


M 

g 


The  main  flow  m'  = (pUA)  = p U sh  CosA,  h. 

g g ^g  g Ch 

2 h 

Coolant  flow  m = 2.48'  . . --  ^c^c 

c X 4 '’l 

Here  the  perimeter  of  tlie  Made  profile  is  taken  as  2.48C,  whicli 
is  equal  to  t lie  perimeter  0!  llie  refi'ience  blade. 

U . V, 

S g throat  X 1 

2.48TI  ■ C ■ t U ■ d ■ d (5.2) 


m 

• = 

lit 


Combining  equations  5.1  and  5 


we  get  the  following  I'xpi'i-ss  ions : 
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For  < 62Z,  i.e..  A,,  < 45 

iso  I 


c , 4.17  m'  -0.35  ^-0.68  , , _-l 

— = 1.944  n — m 0 (s/c  cos  8) 

m m 

8 


where  m'  = 


U 

P c 


PU 


q throat 


(5.3) 


For  33.4%  < q.  ^ < 6-^:  i-e.,  45  A.,  180 

isoth  2 


™c  . 2.273  m'  -0.35  -0.68  ,s 

= 0.833  — ra  (—  cos  6) 

mm  c 

g 


-1 


(5.4) 


For  n . , < 33%,  i.e.,  A„  > 180 

isoth  2 


= 0.380n’"^^^  m'  -0.35  „-0.68  ,s  _ ^ 

m — . m 0 (—  cos  S) 

g m c 


(5.5) 


wliere 


, *(-41  -1.35  ^-0.68  X 

A,  = n*  m 0 


With 


iL-  = A . • Jl 

Seq  7T  d d 


Seg 

for  double  row  of  holes. 


pU 


pU 


m ' _ g throat 

m fiO 

c_ 

(’  . . 
g iniection 

lor  the  relereiu'e  Hl.ide  !. 


((').. 

= injection 

( I'  . 

q throat 


A* 

A.  . 
inject 


0.6  3 


I'.'.ii.it  i('ii  5 1,  v.'hicl)  repri'sents  high  \. lines  of  n.  , , is  used 

' isi'lhermal 

in  this  analysis. 
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Ai-cording  to  the  definition  of  film-cooling  isothermal  effective- 
ness , 


^'isoth 


(5.6) 


where  q 

q' 

q' 


heat  transfer  to  the  blade  in  the  absence  of  film  cooling 

heat  transfer  to  the  blade  with  film  cooling 

(1  - n.  . - and  this  must  be  handled  bv  the  internal 

isothermal 

cooling. 


i ,e.  , 


q"  = m Cp 
c c 


(T 

c . exit 


T 

c 


) 

.m 


(5.7) 


Let  us  define  an  internal  cooling  effectiveness: 


T , - T , 

c exit c . in 

T , - T .in 

blade  c 


(5.8) 


varies  from  0.0  to  1,  higher  values  of  generally  0.5  are  of  interest. 

t is  a strong  function  of  the  type  of  internal  cooling;  its  configuration 
coolant  flow  Reynold  number,  etc. 

Tc  , T - coolant  inlet  and  exit  temperatures  respecitvelv 

in  c exit 

combining  equation  5.7  and  5-8  results  in: 


q"  = in  . C 
^ c pc 


( T - T ) 

^ blade  c.in^ 


^pc 


^ IT.  , , - T . ] 

blade  c . in 


U.  - 'Dll'... 


pc 


[1  - T . 1 

^ bl.idi-  I'.in' 


(1  -■,)  J.  .S  St  . . r . Cp  . K|  (T^ 


T IK 
b 1 


Op 


bl.nli 


wiieft’  St  - blade  mean  Si  ant  on  number 


-43- 


A* 

. (1  - n)  (2.48)a  . St  . CPg  . (T^  - 

c _ exl  t 

“g  Cp  . £ [T,  - T . in]  cos  6 

° b c 


C5.9) 


where  - is  the  coefficient  that  takes  into  account  the  effect  of  blade 
rotation.  Using  experimental  correlation  of  references  (121, 
(13),  and  (14),  is  taken  to  be  1.5. 

2 - exit  angle  at  the  throat. 


m 

1 .mi 

c 

II 

) 

rii 

1 4 ■ mg 

V & 

film  1 

cooling 


internal 
cool ing 


where  K,  - is  the  fraction  of  internal  cooling  being  used  for  film  cooling 
and  taken  to  be  1 in  this  analysis. 

Equations  5.J  and  5.9  are  coupled  and  solved  simul taneousl v . 


= 1.2247n 
mg 


4.17  -0.35  ^-0.68 


(s/c  cos  6) 


-1 


(5.10) 


m (1  - h)  (2.1328)0  . St  . Cp  . K,  [T  - T 1 

c _ ______  S • g 

Cp  . ■ . IT T . ] cos  6 

^c  b c . in 


(^.11) 


i' . ex  i t 

■ Y 


(T  - T . ) + T . 
b c . 1 n c . 1 n 


j :>t  e r n.1  i I I 'o  1 i ng  e i ^ I'l  i i veiu  ss  • . I • = 0 . > , 0 . / ' ,iik1  0.4)  . 
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For  t:ic  pr*,ssure  ratio  transonic  turbine-  stage  tliti 

t -olrnC  ..r.ss  flu---  rc.-quirc-rr.ents  for  tbe  ro'.^r  blav-ie  ro’«  alone  is 

File  coolant  nass  flc-’  rt-qui  at  for  t'ae  norole  row  is  not 
incl-oJeJ  ns  tiie  nczzle  is  con-;  idered  as  part  of  the  ccr.b-astion  chamber. 

The  turbine  inlet  stagnation  absolute  terperature  crfi;>-d  here  is 
for  ti.e  rotor  inlet  Toi,’^  (■=  ibUu°K) 


Blade  Kow : For  the  rotor  blade  row,  equations  5.3  A 5.9  are  solved 


s 1 r.u  1 1 an  c-  ou  s 1 v 
T 


0 = 


c . exi  t 


"b  - h.i„>  + 


,.re  1 


( 5 . ! 2 1 


02 


= 1200  k 


allowable  blade  tcriperature 


"c.in  - 730°k 


...rel  ,--,cO, 

i_-  “ lo75  k 


S/c  cos  B = C0.7)(0.A2)  = 0.315 
o = c/s  =•  ; .4236 
St  - 0.00187 


1.5  accounting  for  the  effect  of  rotation 
„.rel 


E 02 

CoEjiut  u t i ons  aie  i.iiticii  out  tor  r;  =•  0.5,  1.0.  !•  1.5  w;.lie  var>  ing 
( ■ = 0.5,  0.75  and  U.'-) 


a^  Is  defined  as  I h.-  i.stio  of  ronl.ant  to  gar-  flow  rate, 


e r 


; i an  . ii ; - 

St . a 
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Two-Stage  Subsonic  Turbine 

Equation  5.3  and  5.9  are  solved  simultaneously  to  determine  the  coolant 
to  gas  flow  rate. 

As  in  the  transonic  stage;  the  first  nozzle  is  not  included.  This 
is  regarded  as  part  of  the  combustion  chamber. 

The  coolant  mass  flow  requirements  for  the  rotor  of  the  first  stage, 
the  nozzle  and  rotor  of  the  second  stage  are  estimated, 
if  *L  ^ t^ge  Rot ojr 


G is  defined  as  in  equation  5.12 

Tj^  - 1200°K  - allowable  balde  temperature 

T , *■  750°K 

c . in 

T 

02  = T = 1680  K 

S 

(s/c  cos  d)  = (0.7) (0.554)  - 0.387S 

0 •=  c/s  “ 1.4286 

Sj.  = 0.00194 

“ 1.5  accounting  for  the  effect  of  rotation.  Computations  carried 
out  for  n = 0.5,  1.0,  and  1.5  and  t,  “ 0.5,  0.75  and  0.9 
is  defined  as  tiic  ratio  of  coola.nt  to  gas  flow  rate 


1 m 
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2nd  Sta>;e  No zzlc 

The  2nd  stcige  nozzle  inlet  temperature  is  computed  using  energy 
conserva  t ion . 

til  Cp  T‘**^  + ri;  Cp  T 
g c UJ  I . r 1 ' c c exi t 


= (m  -Hit  ,)C'p  'I'^+(m  +1T1  )Cp.AT„m  Cp  + m ,Cp  T 

g c.rl  *q  03  g c.rl  0 g '^g  02  c.rl  i exit 


c . ex  1 1 


= (m  + li,  .'(C’p  ri'''  + C . -T  ) 

g c.rj  g 03  p d 


(5.131 


where  mg  is  the  main  (hot)  gas  mas.s  flow  entering  t lie  first  rotor. 


C.rl 


- coolant  mass  flow  for  rotor  blades  of  Stage  1 


(ih^  + mc.rl)  - mass  flow  leaving  the  rotor  stage 


ATp  - stagnation  temperature  drop  across  the  first  stage. 


l'.t;uat  ion  5.1s  then  reduces  to; 


» ah 
*03 


C P 1 C p . 

'Cp,.  02  1 C'p„ 


T' 


c . ex 1 1 


(1  + a^) 


(5.U) 


ah 

03 


^ 'ii  T-'^  + ^^c.exit  I .... 

C’p  ^02  ^ ‘ ■ ‘c.in] 


(1  + a j ) 


C’p^.  -^01 


knowing 


can  hi'  c.tl ciliated. 


• exit 
I 

k 


IT, 


T . ) + 
I . 1 n 

‘ah 

0 i 


1 , - s t a t ion. 1 1 r ow 


'..'ll  to.  ,S1.’)  (1.4.,,, 
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O = c/s  = 1.4286 

St  = 0.  00194 

ab 

T = T 

g 03  expressed  by  equation  6.2.14 

= 1200“K 
b 

^c.in  = b00°K 

The  equations  are  solved  again  fo.  m = 1.0  and  1.5  and  = 0.5, 

0.75  and  0.9. 

a^  is  defined  as  the  ratio  of  the  coolant  to  gas  flow  rate  for  the 
2nd  stage  nozzle 


m + m 

S ‘^rl 


2nd  Stage  Rotor 

Equations  5.3  and  5.9  are  solved  slroul taneously  to  determine  coolant 
mass  flow. 

„ab  , . . , 

The  mean  absolute  stagnation  temperature  at  2nd  rotor  inlet 

is  found  using  conservation  of  energy. 


(m  + m ,)T,.  , + m ,,  Cp  T . = (m  + m + m (5.17) 

g c.rr  0!  c.n2  c c.exit  g c.rl  c.ii2  4 04 


Cp  , , Cp 

' . T-;';  + a,,  . [T  . 

(.p^  03  ^ c.exit 

(1  + a.,) 


(5. IK) 


'1  .„ab  ^ 

Cp,  ■ ^03  '*2  Cp 

■4  ** 


T ) + T ] 
c . 1 n 1 . 1 n 


(1  + a^) 


re1  _ al. 
Oh  Oh 


2 2 
(C  - K7) 
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where  C'.,  and  W,  - are  the  absolute  and  relative  velocities 
4 4 

(taken  from  the  velocity  triangle) 


0 = 


c . exit 


[,rel 

04 


0 = 


(T  - T . ) + 1 . ] 

b c.in'^  c.in^ 


p3b 

04 


2 

(C,  - W, ) 
4 

2Cp 


i5.20) 


= 1.5 

s/c  Cos  B = (0.7)  x(0.5812)  = 0.4068 

a = c/s  = 1.4286 

St  = 0.00191 

T = compressed  as  in  equation  5.19 

g 04 

T . = 600°K 

c . in 

The  tquations  are  solved  for  m = 0.5,  1.0  and  1.5  and  C = 0.5, 
0.  / ■)  and  0.9. 


m _ , 

a,  is  defined  as  

J m + m , + m 

g c.r  c.n2 


and  a 


subsonii 


. r 1 + m , + , 

c . n 2 i-  . r 2 


ra  t m , + m , + m t 
g c . r I i-  . n2  i-  . r 2 
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The  analysis  yields  the  total  coolant  mass  flows  for  both  the  sub- 
sonic and  transonic  turbines.  The  savings  in  coolant  mass  flow  may 
then  be  computed. 

Limitations  of  the  Anlaysis 

The  major  limitations  of  the  anlysis  lie  with  the  data  used.  For 
example,  the  film-cooling  correlation  used  here  is  one  which  was  obtained 
for  film  cooling  on  a flat  plate.  Thus  the  effects  of  curvature,  the 
effects  of  the  coolant  from  the  preceeding  line  of  holes  and  the  effect 
of  secondary  flows  on  cooling  in  a 3D  machine  are  not  taken  into  account . 
Also  many  other  considerations  such  as  the  losses  due  to  coolant  flow 
injection  and  a detailed  heat  transfer  and  structural  analysis  must  be 
carried  out  before  final  design  criteria  can  be  set. 

5 . 3 Comparison  of  Required  Heat  Removal : 

Without  specifying  any  cooling  technique  to  be  used  a general 
approach  to  determine  the  desirability  of  using  transonic  turbine  to 
replace  subsonic  turbine  can  be  roughly  estimated  by  considering  the 
ratio  of  the  amount  of  heat  that  has  to  be  removed  from  the  stage 
to  the  work  output  from  the  stage. 

5.3.1  Transonic  Stage 

It  is  assumed  that  both  t lu’  transonic  and  subsonic  turbine 
have  the  same  turbine  inlet  total  tenqierature  .ind  total  pressure; 
operate  with  the  same  mass  flow,  axial  velocitx  and  tip  speed.  And 
we  are  considering  Lite  case,  wlu-ri'  the  highlv  lo.ided  transonic  st.ijp' 


Pol  ” 4 is  used  to  replace,  two  stages,  sav  nu'di'ra  t e 1 v Pol 
Po  3 ho  ) 


w 
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subsonic  stages.  For  the  Single  Transonic  Stage, 


transonic 


transonic  n,(C^^  ^ 


where  ^q^^^^g^nic  '^nozzel  + 

blade 

= sum  of  all  the  heat  that  has  to  be  removed  from  the  turbine  components. 
For  simplicity,  let  us  consider  only  the  nozzle  and  blade  rows. 
Using  the  example  we  are  dealing  with  here,  let  us  first  compare 
the  highly  loaded  single  stage  transonic  turbine  with  the  moderately 
loaded  subsonic  one. 

5.3.1  Transonic  Stage : 

For  the  Transonic  Stage  with  a pressure  ratio  of  4,  turbine 
inlet  temperature  of  1800°K  we  have: 


d = 


q , + q 

nozzU'  rotor 


m{;p  AT 


0 


q I q 

nozzle  ^ rotor 


mC  p 


" mC  p 


2.4K  'St  . 1.  i.CpCr  -T  kI 

3..os.st  .K.l.r.tY, i 

1 U (£ 

f 

. U 1 g cos  ■ Cp 

i . r r > • 1 ■ p 

0 
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(2.48) (1.33) (0.000188 )(600) 
0.423 


(1.5)(2.48)(1.333) (0.00188) (1575-1200) 
0.423 


450 


= 0.367 


5.3.2  Subsoni< 


For  the  subsonic  stage,  let  us  consider  the  first  stage  of  the  two- 
stage  turbine. 


d = 


2.48a  St(TQ^-Tj^) 
cos  6 


cos  6 


AT 


01 


= 0.0602 


Second  Subsonic  Stage 

With  pressure  ratio  of  two,  with  turbine  inlet  temperature  of 
15t>0°K 


subsonic  II 


3.7659  + 2.6727 
2'l0° 


= 0.0307 


For  the  Multistage  Subsonic  Turbine 


>1  1,9 

nozzle  ■+■  rotor 


nozz  1 e 


mu  1 1 i St  agt‘ 
subs . 


stage  1 ^ - 

mCp  AT^j  + mCp  AT^„  ^ 


+ rotor 


‘nozzle  I 
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L.e.  by  replacing  two  subsonic  stages  with  a single  highly  loaded 
transonic  stage  to  produce  the  same  work  output,  the  amount  of  licat 
that  has  to  be  removed  from  tlie  turbine  stage  has  been  decreased  by  as 
musch  as  21%. 

It  also  is  wortli  mentioning  that  for  a fixed  rotor  blade  temperature, 
the  turbine  nozzle  inlet  temperature  can  be  raised  for  tlie  transonic  c;ise 
since  tlie  blade  relative  stagnation  temperature  is  generally  lower  in 
transonic  turbine  tlian  in  subsonic  ones. 

i\nd  when  this  increased  temperature  is  matched  with  a <-orresponii  ing 
increase  in  compr  ssor  pressure  ratio,  t'.e  overall  cycle  - etticien>y 
will  improvi’. 

The  analysis  thus  shows  the  potential  lor  gains  in  tin  work  output 
per  unit  ol  looled  blade  surface  tor  liigh  perlormatice  transonic  tut  nines. 


VI.  CONCLUSIONS 


The  most  significant  results  of  the  study  may  be  summarized  in 
the  following  conclusions: 

(a)  Since  the  pressure  distribution  obtained  in  the  cascade  blowdown 
facility  is  the  same  as  that  obtained  in  the  conventional  wind  tunnel 
at  the  Von  Karman  Institute,  then  the  blowdown  facility  is  a bonafide 
and  practical  method  of  aeordynamic  testing.  It  is  also  a flexible 
means  of  solving  rather  complicated  heat-transfer  problems 

while  providing  a fairly  rigorous  modeling  of  the  aerodynamic  flow 
and  temperature  fields  to  simulate  turbine  operating  conditions. 

(b)  A comparison  of  the  isotropic  Mach  number  distribution  obtained 
in  VKI  with  that  obtained  in  the  MIT  cascade  facility  shows  that  the 
level  of  turbulence  in  the  mainstream  has  a marked  effect  on  the 
transition  to  turbulence  in  the  boundary  layer. 

(c)  The  losses  measured  for  transonic  cascades  are  of  the  same  order 
of  magnitude  as  for  subsonic  cascades  which  seems  to  suggest  that 

the  effectiveness  of  both  turbines  are  comparable.  However,  this 
will  have  to  be  modified  by  an  evaluation  of  the  effect  of  coolant 
flow  on  these  efficiencies. 

(d)  A comparison  of  the  loss  curves  for  the  four  blade  profiles  shows 
that  each  profile  has  a superior  performance  in  a different  Mach 
number  range.  Thus  the  reference  blade  with  the  thin  trailing  edge 

is  superior  up  t<'  tiie  design  point  and  tlien  the  plug  nozzle  as  well 
as  the  blade  with  a ciuivergent -d ivergent  nozzle  are  superior  in  t he 
higlier  Mach  number  range.  The  comparison  also  shows  that  the  thickness 
ol  the  trailing  edge  lias  a marked  effect  on  losses. 
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(e)  The  Nusseli  number  distribution  predicted  from  a NASA  computer 
agrees  well  with  the  measurement  except  in  ttie  regions  where  shock 
waves  interact  with  boundary  layers.  This  serves  to  confirm  its  va- 
lidity. The  distribution  of  Nusselt  numbers  also  shows  a very  high 
tieat  transfer  rate  at  the  trailing  edge, 

(f)  The  comparison  of  the  results  obtained  with  those  of  Turner  [15] 
seem  to  suggest  that  the  level  of  turbulence  may  be  a key  variable 

in  determining  heat  transfer. 

(g)  Lastly,  comparative  studies  between  subsonic  and  transonic 
turbines  show  that  there  is  21%  less  heat  to  be  taken  out  by  internal 
cooling  in  the  case  of  transonic  turbines,  and  that  less  coolant 
flow  is  required  to  cool  a single  stage  transonic  turbine  with  a 
pressure  ratio  of  4 used  to  replace  a subsonic  turbine  with  two 
stages  with  pressure  ratios  of  2 per  stage.  All  these  point  to 

the  great  potential  of  transonic  turbines.  But  before  this  potential 
is  realised  further  investigations  are  necessary  such  as  tlic  effect 
of  coolant  flow  Injection  on  the  efficiency  of  transonic  turbines. 


I 
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